Non-polarizing two-atom interaction potential in the liquid 4 He 
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A formula for potential U of the 4 He- 4 He interaction in the liquid state is 
obtained by the direct electromagnetic computation as a function of the in- 
teratomic distance R. The potential decreases exponentially at large R. The 
further development and application of the result are discussed. 
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It is known, that the Leonard- Jones (L-J) potential describes the interaction of atoms 
of simple non-polar gases |]J. However, at large distances between atoms the interaction 
forces fall faster, than it is described by the L-J potential 0. The corresponding cor- 
rections (including light atoms) are suggested in papers 0-[|§. These authors had a 
suspicion, that, for light atoms, either L-J formula or pair interaction approximation, or 
both of them are not valid. In this paper, we show that the rapid decrease of the potential 
can be explained, if we assume that the atom polarization is not important in the liquid. 

In this paper, we calculate the potential U between charge densities of electrons and 
nuclei for 4 He in the ground state. The interaction energy of two charge distributions 
with the densities pi(r) and ^2(1*) is 

U = / , Pi(r)p 2 (r / )drdr / ^ 



where the charge densities are 

pi(r) = q5(r - n) - qP(r - n), 

(2) 

p 2 (r') = qS(r' - r 2 ) - qP(r' - r 2 ). 

Here P(r — r x ) and P(r' — r 2 ) are the electron probability densities at the points r and 
r', the atom nuclei being at the points ri and r 2 , and q = 2e. Using (H) in (|l|), we have 

Pi(r)p 2 (r')drdr' 2 f S(r — ri)S(r' — r 2 )drdr' 

1 ^ 



-q 2 I ^| v 71 I "\ v n + 

I r — r I J I r — r I 

+ ^ P(r - rOPtr' - r 2 )drdr> = ^ + ^ + ^ + ^ 



5(r — ri)P(r' — r 2 )drdr' 2 f 5(r' — r 2 )P(r — r^drdr' 
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It is evident, that the first term U nn in ([!]) gives the interaction energy of two positive 
point-like charges of nuclei after the integration over coordinates r and r': 

q 2 

U = — 
a 

The term U nc describes the interaction between nucleus of an atom and electron cloud 
of another. We have 

_ _ 2 f ^(r-r 1 )P(r / -r 2 )drdr / _ _ 2 f P(r' - r 2 )dr / _ _ 2 r P(r )dr 
J |r — r | J |ri — r | J |K — ro| 

where r = r' — r 2 and R = r x — r 2 . We assume that the electron probability density 
P(i"o) is spherically symmetric, the center being at the nucleus (the accuracy of this 
approximation is discussed at the end of the paper). Integration with respect to the 
spherical coordinates 9 and (p in (HI) gives 



„2 „2 



U nc (R) = + 4vr^ J P(r )(r 2 - Rr )dr (5) 

R 

We used the change of variable 

7r R+ro 

9^y = JR 2 + rl - 2Rr sin9, [ sin 9d9 -> / (6) 
v J J Prn 

|fl-ro| 

The cloud-cloud interaction term [/ cc in (|^) after the similar integration with respect 
to spherical coordinates becomes 

^ R R+x ^ 00 x+R 

U cc = J U nc (x) J P(y)ydyxdx - J U nc (x) J P(y)ydyxdx. (7) 

R-x R x-R 

The formulas for U cc and U nc are valid for any atoms in liquids and gazes at sufficiently 
low temperatures. Concrete expressions for P(r) have to be used in (U) and (0) to obtain 
the final formula. In particular, for 4 He, both electrons are at the state with n — 1, 1 = 
(para-state). To obtain a simple formula for using in preliminary computations, we may 
ignore the electron-electron interactions in the cloud. Then the wave function of the 
electrons in the 4 He ground state is 

$(r 3 ,r 4 ) = ^(r 3 )^(r 4 )x P , 
where r 3 , r 4 are the electron coordinates, 

* (r) = ^r r/ °' * = >T«-UT>], 

X p is antisymmetric spin wave function of the para-state. 

The product ^/(r 3 )^/(r4) describes the symmetric radial part of the ground-state wave 
function. |$(r 3 ,r4)| 2 is the probability density to find one electron at the point r 3 , and 
another at the point r 4 . Because \x P \ 2 = 1, the probability density is 

|$(r 3 ,r 4 )| 2 = ^ 2 (r 3 )^ 2 (r 4 ). 
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The charge density created by the two-electron cloud at the point r is 



qP(r) = e|^(r)| 2 + e|^(r)| 2 = g|*(r)| 2 . 
Using this formula and (|J), (0), we obtain the expressions for U nc and U cc : 

1 V 



R 



U n c ' 
1 1 



2R/a 



V/ 



a. 



6a 



Thus, the interaction potential between two atoms at a distance R is 

q 2 \ 1 /R\ 2 3 fR\ 5 



U(R) = U nn + + U nc + U CC = — 

a 



6 V a ) 



4 \a) + 8 + Vi? 



-2-R/a 



(9) 



(10) 



It is obvious, that 



U(R) 



q 2 R 2 



- 2R/a , as R > a. 



The exponential decrease of our potential U contrasts with i?~ 6 -decrease in other formulas 
for 4 He- 4 He potential. The slow i? -6 -decrease is provided by the use only the dipole 
interaction is taken into account in these formulas. 

Our formula ( |I~0"D for U (R) does not contain arbitrary parameters. It is approximate, 
however, because we neglected the electron-electron interaction for obtaining ^(r^,r4). 
Nevertheless, formula suggest an idea, that the semi-phenomenological potential of the 
form 



U(R) 



--yR/a 



(11) 



is of interest for the theory of the liquid 4 He maxon-roton spectrum, that considerably 
depends on the form of interatomic potential 0. 

To take the electron interaction in the cloud into account, we have to use some realistic 
wave functions of the 4 He electrons in the ground state (see, e. g., M). For instance, the 
wave function of Hylleraas with three parameters gives additional terms up to i? 6 inside 
the square brackets in (|TT|) . We emphasize that our formulas are designed for describing 
the two-atom interactions in the liquids and gases and we suppose, that no polarization 
effect exists for any pair of atoms in liquid due to its compensation. 

The interaction potential of two atom in the vacuum in general has a contribution 
caused by the atom polarization. The atom polarization results in that the centers of 
electron clouds do not coincide with the nuclei and the clouds are deformed. The first 
estimate of a polarization effect can be easily described by our formulas with R in U nc 
less than R in U nn and U cc - In such case, taking A as a shift of a positive and negative 
charge centers, we have a Buckingham-type formula 



U(R) = Unn(R) + 2U nc (R ~ A) + U CC {R) 



2q 2 A 



R(R-A) 



+ V(R), 



(12) 



where V(R) decreases as exp(— 2R/a) as R 3> a. It is important, that A depends on 
R and A — ► as (R/a) — > oo. In other words, the dipole moment caused by the atom 
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polarization decreases as (R/a) — > oo. It is very difficult to specify the dipole moment 
dependence on R. Nevertheless, we can say that i? -6 -term in the L-J-potential is not 
correct in liquid media. As to the pair of 4 He atoms in vacuum, R~ 6 term is not correct 
due to the dependence A(R) (that is unknown unfortunately), and all of these properties 
are inherent for 4 He as the most compact inert gas. 
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